A new wall-damping function for large eddy simulation employing Kolmogorov velocity scale

A new wall-damping function, based on the Kolmogorov velocity scale, for large eddy simulation (LES) is proposed, which accounts for the near-wall effect. To calculate the Kolmogorov velocity scale, u_ε, the dissipation rate of turbulent energy, ε, is needed. In LES, however, the dissipation rate is generally not solved, unlike in the Reynolds averaged Navier-Stokes (RANS) simulations, e.g., k-ε models. Although, in some previous studies, the dissipation rate of the subgrid-scale (SGS) turbulent energy, ε_SGS, is used instead of ε in calculating the Kolmogorov velocity scale, the scale obtained using such a method overly depends on the grid resolution employed and is generally inappropriate. Accordingly, the wall-damping function using the incorrect velocity scale also depends on the grid resolution and gives an inadequate wall effect. This is because ε_SGS contains only the components in the scale smaller than the grid-filter width, which obviously varies with the grid resolution employed. In this study, to overcome this problem, we propose a method for estimating the Kolmogorov velocity scale with a technique of conversion in LES, and the estimated one is utilized in the wall-damping function. The revised wall-damping function for LES is tested in channel flows and a backward-facing step flow. The results show that it yields a proper near-wall effect in all test cases which cover a wide range of grid resolution and Reynolds numbers. It is also shown that all three kinds of SGS models incorporating the present wall-damping function provide good predictions, and it is effective both in one-equation and 0-equation SGS models. These results suggest that the use of the proposed wall-damping function is a refined and versatile near-wall treatment in LES with various kinds of SGS models.     

Effect of thermophoresis and other parameters on the particle deposition on a tilted surface

In this study a modified version of v2-f turbulence model (φ-α), is applied to simulate a non-isothermal air-flow. The φ-α model and a two-phase Eulerian approach complement each other to predict the rate of particle deposition on a tilted surface. The φ-α model can accurately calculate the normal fluctuations, which mainly represent the non-isotropic nature of turbulence regime near the wall. The Eulerian model was modified considering the most important mechanism in the particle deposition rate when compared to the experimental data. The model performance is examined by comparing the rate of particle deposition on a vertical surface with the experimental data in a turbulent channel flow available in the literature. The effects of lift force, turbophoretic force, thermophoreric force, electrostatic force, gravitational force and Brownian/turbulent diffusion were examined on the particle deposition rate. The results show that, using the φ-α model predicts the rate of deposition with reasonable accuracy. The results of modified particle model are in good agreement with the experimental data. This study highlights the paramount effect of thermophoretic force on the particle deposition rate and clearly shows that when the temperature difference exceeds a certain limit, the electrostatic force has insignificant effect on the particle deposition rate. Furthermore, it is indicated that even at small temperature differences, the effect of tilt angle on the particle deposition rate for intermediate-size particles is negligible.     

The issues of the uniqueness and the stability of the homogeneous response in uniaxial tests with gradient damage models

We consider a wide class of gradient damage models which are characterized by two constitutive functions after a normalization of the scalar damage parameter. The evolution problem is formulated following a variational approach based on the principles of irreversibility, stability and energy balance. Applied to a monotonically increasing traction test of a one-dimensional bar, we consider the homogeneous response where both the strain and the damage fields are uniform in space. In the case of a softening behavior, we show that the homogeneous state of the bar at a given time is stable provided that the length of the bar is less than a state dependent critical value and unstable otherwise. However, we also show that bifurcations can appear even if the homogeneous state is stable. All these results are obtained in a closed form. Finally, we propose a practical method to identify the two constitutive functions. This method is based on the measure of the homogeneous response in a situation where this response is stable without possibility of bifurcation, and on a procedure which gives the opportunity to detect its loss of stability. All the theoretical analyses are illustrated by examples.     

On the predictive capabilities of the shear modified Gurson and the modified Mohr-Coulomb fracture models over a wide range of stress triaxialities and Lode angles

The predictive capabilities of the shear-modified Gurson model [Nielsen and Tvergaard, Eng. Fract. Mech. 77, 2010] and the Modified Mohr-Coulomb (MMC) fracture model [Bai and Wierzbicki, Int. J. Fract. 161, 2010] are evaluated. Both phenomenological fracture models are physics-inspired and take the effect of the first and third stress tensor invariants into account in predicting the onset of ductile fracture. The MMC model is based on the assumption that the initiation of fracture is determined by a critical stress state, while the shear-modified Gurson model assumes void growth as the governing mechanism. Fracture experiments on TRIP-assisted steel sheets covering a wide range of stress states (from shear to equibiaxial tension) are used to calibrate and validate these models. The model accuracy is quantified based on the predictions of the displacement to fracture for experiments which have not been used for calibration. It is found that the MMC model predictions agree well with all experiments (less than 4% error), while less accurate predictions are observed for the shear-modified Gurson model. A comparison of plots of the strain to fracture as a function of the stress triaxiality and the normalized third invariant reveals significant differences between the two models except within the vicinity of stress states that have been used for calibration.     

Micromechanical modelling of remanent properties of morphotropic PZT

The present paper investigates the capability of micromechanical material models to predict the ferroelectric behaviour of morphotropic PZT ceramics in a rate-independent approximation based on realistic microscopic material parameters. Starting point is a three-dimensional tetragonal model, which builds on the model of Pathak and McMeeking [2008. Three-dimensional finite element simulations of ferroelectric polycrystals under electrical and mechanical loading. Journal of the Mechanics and Physics of Solids 56, 663-683]. Volume fractions of the crystallographic variants represent the domain structure inside the grains. Interactions between the grains are taken into account by means of a representative volume element of the grain compound. A simplified set of realistic microscopic material parameters of the lattice in terms of Young's modulus, Poisson's ratio, dielectric constant, and spontaneous strain and polarisation is derived from experimental data and theoretical results given in the literature. The simulation of the macroscopic remanent polarisation and strain response due to two load cases shows explicitly that the tetragonal model is not capable to reproduce the behaviour of morphotropic PZT. Therefore, the model is extended by the rhombohedral phase, allowing a mixture of both phases with varying quantities inside the grains. A comparison of our results with experimental data shows a remarkably good agreement, revealing the capability of the extended model.     

The interlayer shear effect on graphene multilayer resonators

Graphene nanostrips with single or few layers can be used as bending resonators with extremely high sensitivity to environmental changes. In this paper we report molecular dynamics (MD) simulation results on the fundamental and secondary resonant frequencies f of cantilever graphene nanostrips with different layer number n and different nanostrip length L. The results deviate significantly from the prediction of not only the Euler-Bernoulli beam theory (f∝nL⁻²), but also the Timoshenko's model. Since graphene nanostrips have extremely high intralayer Young's modulus and ultralow interlayer shear modulus, we propose a multibeam shear model (MBSM) that neglects the intralayer stretch but accounts for the interlayer shear. The MBSM prediction of the fundamental and secondary resonant frequencies f can be well expressed in the form f−f_mono∝[(n-1)/n]^bL^−2(1−b), where f_mono denotes the corresponding resonant frequency as the layer number is 1, with b=0.61 and 0.77 for the fundamental and secondary resonant modes. Without any additional parameters fitting, the prediction from MBSM agrees excellently with the MD simulation results. The model is thus of importance for designing multilayer graphene nanostrips based applications, such as resonators, sensors and actuators, where interlayer shear has apparent impacts on the mechanical deformation, vibration and energy dissipation processes therein.     

侧伏规律在岩质边坡结构面编录中的应用

理论上,据边坡设计坡面来编录展示结构面的迹线分布应与其在坡面内的侧伏交线一致。把握一定的侧伏规律,处以适当的素描技法,既可吻合实际情况,又可提高绘图效率和资料质量,对边坡结构模型的准确分析不无裨益。有鉴于此,经过现场观测总结,分析了结构面与坡面的各种交切关系及其迹线的侧伏特征。在侧伏规律公式推导的基础上,根据两者产状提出其特征三角形判别法,并编制相应的算图,可校核以及在适当条件下推测某一产状要素。讨论了实践中几类难以作图的特殊情况下侧伏规律的应用。     

On stress-state dependent plasticity modeling: Significance of the hydrostatic stress, the third invariant of stress deviator and the non-associated flow rule

It has been shown that the plastic response of many materials, including some metallic alloys, depends on the stress state. In this paper, we describe a plasticity model for isotropic materials, which is a function of the hydrostatic stress as well as the second and third invariants of the stress deviator, and present its finite element implementation, including integration of the constitutive equations using the backward Euler method and formulation of the consistent tangent moduli. Special attention is paid for the adoption of the non-associated flow rule. As an application, this model is calibrated and verified for a 5083 aluminum alloy. Furthermore, the Gurson-Tvergaard-Needleman porous plasticity model, which is widely used to simulate the void growth process of ductile fracture, is extended to include the effects of hydrostatic stress and the third invariant of stress deviator on the matrix material.     

Dual-stage nested homogenization for rate-dependent anisotropic elasto-plasticity model of dendritic cast aluminum alloys

This paper proposes a nested dual-stage homogenization method for developing microstructure based continuum elasto-viscoplastic models for large secondary dendrite arm spacing or SDAS cast aluminum alloys. Microstructures of these alloys are characterized by extremely inhomogeneous distribution of inclusions along the dendrite cell boundaries. Traditional single-step homogenization methods are not suitable for this type of microstructure due to the size of the representative volume element (RVE) and the associated computations required for micromechanical analyses. To circumvent this limitation, two distinct RVE’s or statistically equivalent RVE’s are identified, corresponding to the inherent scales of inhomogeneity in the microstructure. The homogenization is performed in multiple stages for each of the RVE’s identified. The macroscopic behavior is described by a rate-dependent, anisotropic homogenization based continuum plasticity (HCP) model. Anisotropy and viscoplastic parameters in the HCP model are calibrated from homogenization of micro-variables for the different RVE’s. These parameters are dependent on microstructural features such as morphology and distribution of different phases. The uniqueness of the nested two-stage homogenization is that it enables evaluation of the overall homogenized model parameters of the cast alloy from limited experimental data, but also material parameters of constituents like inter-dendritic phase and pure aluminum matrix. The capabilities of the HCP model are demonstrated for a cast aluminum alloy AS7GU having a SDAS of 30 μm.